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Ensemble in-equivalence in supernova matter within a simple model 
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A simple, exactly solvable statistical model is presented for the description of baryonic matter in 
the thermodynamic conditions associated to the evolution of core-collapsing supernova. It is shown 
that the model presents a first order phase transition in the grandcanonical ensemble which is not 
observed in the canonical ensemble. Similar to other model systems studied in condensed matter 
physics, this ensemble in-equivalence is accompanied by negative susceptibility and discontinuities 
in the intensive observables conjugated to the order parameter. This peculiar behavior originates 
from the fact that baryonic matter is subject to attractive short range strong forces as well as 
repulsive long range electromagnetic interactions, partially screened by a background of electrons. 
As such, it is expected in any theoretical treatment of nuclear matter in the stellar environement. 
Consequences for the phenomenology of supernova dynamics are drawn. 

PACS numbers: 64.10. -|-h, 64.60.-i, 26.50.-l-x, 26.60.-c 



I. INTRODUCTION 



Standard thermodynamics is based on the assumption that the physical properties of a system at equilibrium do not 
depend on the statistical ensemble which is used to describe it. Under this condition thermodynamics is unique and 
^ the different thermodynamic potentials are related via simple linear Legendre transforms. Is is however well known 
(-H I that the equivalence between the different statistical ensembles can only be proved [l[ at the thermodynamic limit 
and under the hypothesis of short range interactions, while non-standard thermostatistic tools have been developed 
during the years to deal with non-extensive and long-range interacting systems 0, Q . 
^ I The issue of ensemble in-equivalence, namely the possible dependence of the observed physics on the externally 
' applied constraints, has been typically associated to phase transitions, more precisely to phase separation quenching 
CO . due to the external constraint. A well-known example in the literature concerns the possible occurrence of negative 
' heat capacity in finite systems, which has been widely studied theoretically Q and has also given rise to different 
. experimental applications in nuclear and cluster physics [Hj . In this specific example the phase separation is quenched 
' by the microcanonical conservation constraints, leading to the thermodynamic anomaly of a non-monotonous equation 
of state. 

Concerning macroscopic systems, different model applications have shown fingerprints of ensemble in-equivalence 
but phenomenological applications arc scarce. In this paper we show that the dense matter which is produced 
in the explosion of core-collapse supernova and in neutron stars is an example of a physical system which displays 
this in-equi valence. 

5-H ' We will limit our discussion to finite temperature T « 10^*^ A' and nuclear sub-saturation 10^*^ < p < 10^'* g cm^^ 
densities, thermodynamic conditions which are known to be largely explored in the dynamics of supernova matter 



and in the cooling phase of proto- neutron stars [111 Il2| . The baryonic component of this stellar matter is given by 
a statistical equilibrium of neutrons and protons, the electric charge of the latter being screened by an homogeneous 
electron background. 

If the electromagnetic interactions are ignored, this gives the standard model of nuclear matter, which is known 
to exhibit first and second order phase transitions with baryonic density as an order parameter, meaning that the 
transition concerns a separation between a dense (ordered) and a diluted (disordered) phase [l3|. It is however 
known since decades to the astrophysical community that the situation is drastically different in stellar matter, where 
microscopic dishomogeneities are predicted at almost all values of temperature, density and proton fraction and 



thermodynamical quantities continuously change at the phase transition IJ, [15|. This specific situation of stellar 
matter respect to ordinary nuclear matter has been shown to be due to Coulomb frustration which quenches the first 
order phase transition p^j. However, the thermodynamic consequences of this specific thermodynamics with long 
range interactions have never been addressed to our knowledge. 

In this work we will show that these dishomogeneities imply ensemble in-equivalence, making neutron star matter 
the first astrophysical example to our knowledge of ensemble in-equivalence at the thermodynamic limit. Additionally, 
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we will show that a consistent treatment of this specific thermodynamics can have sizeable effects in the equations of 
state which are currently used to describe the supernova phenomenology. 

In a recent paper [l7j we have proposed a phenomenological hybrid model for supernova matter which was numeri- 
cally solved by Monte-Carlo simulations. Since convergence is always an issue in Monte-Carlo calculations, we propose 
in this paper an analytic version of the same model. In order to have analytical results, we will limit to the simplest 
version of the model where a schematic nuclear energy functional is used. It is clear that more sophisticated energy 
functionals will have to be implemented in order to have quantitative predictions for the supernova simulations. 

However, the qualitative conclusions of this paper will only depend on the sign (attractive or repulsive) of the 
interactions and on the size dependence (in terms of volume and surface) of the nuclear binding energy. As such, they 
will not depend on the details of the model. 



II. THE MODEL 

Stellar matter at temperature lower than the typical nuclear binding energy Ei, ^ 8 MeV/nucleon « 10^^ K, and 
density p lower than the saturation density of nuclear matter po w 0.16 fm~'^ w 1.6 • 10^^ g-cm~^, can be viewed 
as a statistical mixture of free protons and neutrons with loosely interacting nuclear clusters at internal density po, 
immersed in an homogeneous electron background density pe which neutralizes the total positive charge density pp 
over macroscopic length scales, pe = Pp- Nucleons bound in clusters can be described by a phenomenological free 
energy functional depending on the cluster size a = n + z and chemical composition i = n — 2 as well as on the 
temperature of the medium: 

fa,^ = ea,^ + {e^d P ~ T S^^^ (1) 

For nuclear (fermionic) clusters both the average cluster excitation energy (e* and entropy ^ can be evaluated 
in the low temperature Fermi gas approximation 

(e;,)^ = coaT^; (2) 
sf,, = (2coT + csa-^/3/i(r))a, (3) 

The surface term csa^^'^h{T) in Eq. (jSj effectively accounts for the entropy increase at finite temperature due to 
surface excitations, producing a vanishing surface free energy at a given temperature, corresponding to the critical 
point Pc = of nuclear matter. The cluster energy ea,i is modified with respect to the energy of the nucleus in 
the vacuum ^ because of the electromagnetic interaction with the electron background which neutralizes the proton 
charge. In the Wigner-Seitz approximation the functional results [13] 



where pp = pe is the proton density and pop ^ Pp the corresponding saturation value. The minimal scale at which 
charge neutrality is verified, is called a Wigner-Seitz cell. 

We use in the following a simple liquid-drop parameterization for the cluster energies e° ^ [l7j : 

6°,^ = {~cva + csa2/3^ j^i _ ^/^^ + ccz^a-^/^ (5) 

but shell and pairing corrections can be readily incorporated figj . Density dependent correction terms accounting for 
the nuclear interaction with the free nucleons are also expected to improve the predictive power of the model. 



as well as a more sophisticated form for the cluster internal entropy using realistic a and i dependent densities of 
states [TtI- Since none of these improvements is expected to change the qualitative results of this paper, we are not 
including them here to keep an analytically solvable model. 

As a first approximation, one can consider that the the system of interacting nucleons is equivalent to a system 
of non- interacting clusters, nuclear interaction being completely exhausted by clustcrizationp^ . This classical model 
of clustcrized nuclear matter is known in the literature as nuclear statistical equilibrium (NSE) (25l - [28| . This simple 
model can only describe diluted matter at p <^ po a-s it can be found in the outer crust of neutron stars, while nuclear 
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interaction among nuclcons and clusters has to be included for applications at higher density, when the average 
inter-particle distance becomes comparable to the range of the force. 

In our model, interactions among composite clusters are taken into account in the simplified form of a hard sphere 
excluded volume. Since the nuclear density is (approximately) constant inside the clusters, the volume occupied by 
each species {a,i) is given by Va^i = aua^i/po, where Ua^i is the abundance of the species (a, i). The volume fraction 
available to the clusters then reads: 



Vf Pel 
— = 1 — > a — — = 1 

where pd is the total density of nucleons bound in clusters. A slightly different expression has been used in Refs 



V poV po ' 



171 [19|, which however does not change the results presented in this paper. 

In addition to the excluded volume effect for the clusters, the inter-particle nuclear interaction is also considered 
for the nucleons not bound in clusters. The free nucleons self-energy is computed in the self-consistent Hartree-Fock 
approximation with a phenomenological realistic effective interaction [29| . The energy density can be expressed as a 
function of the density of neutrons p„ and protons pp which are not bound in clusters as : 

-f to{xQ + 2)(p„ + ppf/A - tn{2xo + l)ipl + pI)/4: 

+ t3{x3 + 2){pn + PpY+^/2A - t3(2a-3 + l){pn + PpV {pl + pI)/2A 

+ {ti{xi+2)+t2{x2+2)){pn+ Pp){Tn+Tp)/d>+{t2{2x2 + l)-ti{2xi + l)){pnTn+ PpTp)/8, (7) 

where mo is the nucleon mass, to,ti,t2,t3,xo,xi,X2-LX3,a are Skyrme parameters and T„,Tp represent the neutron 
and proton kinetic energy density. We use the SKM* [29| parameterization for the numerical applications. 

At density p > po the whole baryonic matter is expected to be homogeneous and described by eq. ^ . The transition 
from inhomogeneous to homogeneous matter is physically realized in stellar matter at the interface between the crust 
and the core of a neutron star. As we can see from eq.©, in our model at proton densities Pp = pop the Wigner- 
Seitz correction exactly compensates the Coulomb self-energy of the cluster and the total Coulomb energy vanishes, 
reflecting matter homogeneity at supersaturation densities. In this regime the asymptotic cluster energy represents 
the energy density of homogeneous neutral nuclear matter 

lim -ea,i{pp = Pop) = {Pn,Pop), (8) 

a-i-oo a p 

where p = limA,v-i.oo A/V = Pn + Pp, Pi = liui/,\/_>.oo I /V = Pn ~ Pp-, and A, I are the total number of particles and 
chemical asymmetry. The crust-core transition can thus be seen equivalently as the melting of clusters inside dense 
homogeneous matter, or as the emergence of a percolating cluster of infinite size. 

Charge neutrality is imposed globally but local charge dishomogeneities at the scale of the Wigner-Seitz cell naturally 
appear in the thermodynamic conditions where matter is clusterized. In turn, this gives rise to a long range monopole 
component of the Coulomb potential which extends over domains of the order of the cluster size, and which can 
potentially become macroscopic in the limit of very extended clusters close to the crust-core transition. As we will 
show in detail, these long range Coulomb correlations are at the origin of the specific thermodynamics. 



Grandcanonical formulation 



Considering that the center of mass of composite fragments can be treated as a classical degree of freedom, their 
grandcanonical partition sum reads [H [H 



2'a>l 



W exp Za-. 

a> l,2e(— a,a) 



(9) 



The partition sum associated to a cluster composed of n neutrons and z protons in a volume V is given by 
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exp 



(10) 
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Here, Vp is the free volume associated to the cluster center of mass given by eq.®, and the cluster chemical potential 
is a linear combination of the isoscalar and isovector chemical potentials ^, fij which have to be introduced in the 
presence of two conserved charges 



(11) 



Fermi statistics cannot be neglected when treating a = 1 fragments (protons and neutrons). This component of the 
baryonic partition sum is thus included in the finite temperature Hartree-Fock approximation 17 [ 
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where the non-interacting part of the partition sum can be expressed as a functional of the kinetic energy density 
for neutrons {q = n) or protons (q = p) : 



2/3V 



q—n,p 



dTn 



(13) 



with: 



8. 1 p^dp ^^^^ 



' Jo ft2 1 + e^('=.-^.) ■ 

The number density of free protons and neutrons pfq = dp^i^^'^Zp^^^jV determines the single particle energy 
which enters in the self-consistency equation (|14p according to, 

Q^HM Q^HM 

eq=P^^ + . 15 

9rfq Pfq 

The computation of all thermodynamical variables is straightforward using the standard grandcanonical expressions 
for the global baryonic partition sum Z ~ Z°'^^Z°'^^. In particular the total baryonic pressure is simply /3p = 
XvlZp ^^ ^^JY, the multiplicity of the different clusters is given by 

din Za n ai 

ria.i = r^n = Za,i (16) 

and the total baryonic density is the sum of the clusterized (a > 1) and the unbound (a = 1) component 

Pp + Pn = Y^ + P! = Pel + Pi (17) 

where pf = pfp + pf^ is the total density associated to the nucleons which are not bound in clusters. 

Expression Eq. Q includes the electrons self-energy inside the Wigner-Seitz cell. This means that in the global 
stellar partition sum the remaining electron contribution will be factorized Ztot = ^jTJj. ' where Z^'^J^ is given 

by Eq. ([9|) and Z'p is a trivial ideal Fermi gas contribution which has no influence on the thermodynamics and will 
not be discussed further flij . 

The functional relation between the different grandcanonical variables is represented in Fig. [l]for a representative 
given value of T and pi, relevant for the astrophysical applications. The constrained entropy per baryon s^g.^^ = 
'^p.p.il P is evaluated from the Legendre transform of the partition sum in the region where the density is defined: 

(p) = In Zft^^,^^,jV - Ppp, (18) 

and shown in panel (a) as a function of the volume per baryon or inverse density v = 1/ p. The first derivative of this 
function is the baryonic pressure P given in panel (b). The first derivative of the entropy density ap^p,j with respect 
to the density gives the other equation of state, namely the chemical potential p = P/ P — sp^^^/fi, represented in 
panel (c). Finally phase equilibrium is best spotted by looking at the phase diagram given by the relation between 
intensive variables, as shown in panel (d). 
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FIG. 1: Grandcanonical thermodynamics at T=1.6 MeV and /^/=1.68 MeV. Constrained entropy per baryon (a) and associated 
pressure (b) as a function of the inverse baryonic density; to better evidence the concave behavior of s^^^^- a hnear function 
{I3Pt/p — PfJ-ir) is subtracted; (c) chemical potential as a function of the baryonic density; (d) pressure as a function of the 
chemical potential. Dashed lines: Gibbs construction. 



The grandcanonical thermodynamics leads to a first order phase transition. This can be inferred from the charac- 
teristic two-humped structure of the constrained entropy, which is better evidenced by subtraction of a straight line, 
as well as by the crossing of the two equations of state (panel (d)). In the presence of a first order phase transition 
only in the ensemble where the order parameter is fixed (here: the canonical ensemble) the phase transition region 
is accessible, while it is jumped over in the ensemble where the order parameter is fixed only in average (here: the 
grandcanonical). As a consequence, a discontinuity is observed in the equations of state covering a huge range of 
baryonic densities relevant for the description of supernova dynamics. 

The equilibrium solution in the phase transition region corresponds to a linear combination of the two pure phases 
following Gibbs rules. This construction is exactly equivalent to a one-dimensional Maxwell construction if we work 
in an ensemble where all intensive parameters are fixed but one, and it is represented by the dashed lines in Fig. 
[1] (sol . l3l| . We note on passing that this simplification demands to work in the non-standard ensemble (T, p,^/). 
In astrophysical applications it is customary to work rather with the parameters (T,p,yp), where = Pp/p is the 
proton fraction. Within this ensemble, a full two-dimensional Gibbs construction is needed to correctly calculate the 
coexistence zone. The use of a Maxwell construction in the ensemble (T, p, y^) is a mistake, yet it is often used in the 
literature 
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The construction of a convex entropy envelope allows to recognize that inside the density region, indicated by dotted 
lines, the obtained grandcanonical solutions do not correspond to an equilibrium, since a higher entropy solution can 
be obtained by making a linear combination of the two (dense and diluted) solutions marked by a filled circle, that is 
by a first order phase transition. The pressure and chemical potential associated to the metastable branch on the low 
density side are also represented in Fig. [T] From the pressure viewpoint (panel (b)) such solutions appear equivalent 
to a phase coexistence (dashed line), but this is not true in terms of chemical potential (panel (c)). 



The presence of metastable solutions was never discussed in the framework of NSE models [19|, |25| - |28| | to our 
knowledge. It can be understood from the fact that, for chemical potentials higher than the Fermi energy of dense 
uniform matter fi > ep ~ —IGMeV, the equilibrium condition can be obtained either as a mixture of clusters and 
homogeneously distributed nucleons ^"^^ = = fi, or alternatively letting the clusterized component to vanish 

(^°=i = ^ and Pel = 0). This gives a second stationary entropy solution which appears to correspond to the absolute 
entropy maximum, and which renders metastable the clusterized solution at lower density. 

The ending point of the metastable branch can be worked out easily. Indeed from Eq. (jlOp we can see that, for any 
temperature /3~^ and isovector chemical potential /x/, it exist a limiting value of chemical potential ^max for which 
the cluster multiplicities uaj asymptotically diverge. This value is determined by the condition of an asymptotically 
negative Gibbs free energy 



lim /f , -fla-^lI^<0 (19) 



combined with the requirement that the Coulomb energy vanishes in the limit of an infinitely extended homogeneous 
cluster: 

1 \ ^ a — i , 

-riaA = Pop (20) 



V ^ 2 



For chemical potentials above ^imax, clusterized partitions become unstable due to the emergence of a liquid phase. 

Coming back to the discontinuity shown in the equations of state in Figure (TJ in standard thermodynamics the fact 
that the Gibbs construction is a posteriori made to fill up the coexistence region is not a limitation, because ensemble 
equivalence guarantees that the very same linear combination solution would have been obtained if we had worked in 
the canonical ensemble, explicitly constraining the density. 

Because of the complexity of the phenomenology of stellar matter we have not proposed an Hamiltonian treatment of 
the problem. However the phenomenological free energy functional Eq. ([5]) we have used implicitly contains the effect 
of the attractive short range nuclear forces, scaling proportionally to the number of particles in the thermodynamic 
limit, and repulsive long range Coulombic forces, scaling proportionally to the square of the number of particles and 
only partially screened. None of the possible improvements on the nuclear energy functional would change these very 
general scaling behaviors. It has been argued in ref. that this generic frustration effect should lead to a quenching of 
the first order phase transition. We have shown that the phase transition is observed in the grandcanonical ensemble. 
We turn therefore to explore the possibility that ensemble in-equivalence might be at play in the stellar environement, 
with the thermodynamic anomalies associated to the thermodynamics of long range interactions 0] . 



Canonical formulation 



A fully canonical formulation of our model would imply the use of two independent extensive variables, the proton 
and neutron density {pp,pn) or equivalently the isoscalar and isovector density {p,pi). It is however well established 



3^|33| that, contrary to other physical systems like binary alloys and molecular mixtures where phase transitions can 
also imply separation of the species [s^ ] , the isovector density is not an order parameter of the nuclear matter phase 
transition. Since the in-equivalence effect we arc looking for is associated to the phenomenon of phase coexistence, 
we thus expect to see it even if we keep a grandcanonical treatment for the isovector density. 

The effect we are looking for is due to the long range Coulomb interaction, which vanishes in homogeneous matter. 
This Coulomb interaction is neglected in the computation of the abundances of a = 1 particles which are modelized 
as homogeneously distributed. For this reason we will also stick to a grandcanonical formulation for a = 1 particles 



according to Eq. ([T2|) , and assume that the approximate Legendre transformation Eq. ([T8| 

In (p) = In Z^=i^^ - P^xpV, (21) 

which is exact in the mean- field approximation we have employed, is physically correct. For the same reason we will 
consider that the standard thermodynamic assumption of total canonical entropy being additive among independent 
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components, is verified as expected at the thermodynamic limit within ensemble equivalence: 

a^^'^-^ (p) = In Z°,-\ {pf ) + ^lini^ 1 In Z^^>] {Vp,i) (22) 

where the density repartition between the clustered pd and unbound pf component Eq. (|17p is uniquely defined by 
the condition of having a single chemical potential p for both components. Possible extra deviations from ensemble 
equivalence originating from the a = 1 contribution would need a more sophisticated model where the polarization 
of the free protons would be explicitly accounted for [l^ . For simplicity, in the following we will refer to this hybrid 
{P,pi,p) ensemble as to the "canonical" ensemble. 

To derive the expression of the canonical partition sum, we start from the general statistical mechanics relation 
which links the different statistical ensembles, restricted to composite clusters a > 1 only (the subscript a > 1 is 
omitted hereafter for simplicity): 

Zf!,^,^, = ^/3./.. (^) cxp f3pA (23) 

A>1 



Identification with Eq. ([9|) gives. 



%..(^)=En7^ (24) 

{n„}a>l 



where rij, = Uia is the occupation number of size a and the sum is restricted to combinations {ua} = {ri2, ■ . ■ , 'T-a} 
satisfying the canonical constraint, 



A 



CiUa = A. (25) 

The weight of the different cluster size is given by, 

a 

= 0Ja,a-2z exp{l3pi{a - 2z)), (26) 



OJa 

z=0 



and the weight of each nuclear species is determined by the free energy functional we have assumed. 

Following Ref. (35| we introduce an auxiliary canonical partition sum -Z™^^ (A) defined by the additional constraint 
that the cluster multiplicity is fixed to m, J^a "-a ~ ™- With the help of relation we get, 

^;:;/(^ - «) = (28) 

where < Ua >m is the average multiplicity of size a under the additional constraint of total multiplicity m. Explicitly 
implementing the canonical constraint Eq. (j25p we arrive to the recursion relation. 



a=2 

which is valid for any value of A. Summing over all the possible m values, we can see that the same relation holds 
for the unconstrained canonical partition sum 

1 

(^) = ^ E a^aZp^^, {A - a). (30) 

a=2 

This expression can be recursively computed with the initial condition Zi3^^j{2) = a;2- The computational problem 
however arises that going towards the thermodynamic limit the evaluation of the double sum implied by Eq. (|26p be- 
comes numerically very heavy. Starting from a sufficiently high minimal value a„iin we therefore develop a continuous 
approximation for Eq. (j26|) . 

1 r 

Wa>a„.„ ~ 2 / dycyi^g{y), (31) 
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This integral is calculated in the saddle point approximation: 



= lnwa,i + Pi^ii « lni:ja,<i> + Pui < i > r-j ■ (32) 



where the most probable isotopic composition < i > of a cluster of size a depends on the temperature according to, 



and the associated dispersion is given by 



The global weight of size a finally results, 



dfa i 

IJ-i = \i=<%>i (33) 



v/27rcr2exp(/3/^/ <i>)/2 (35) 



is chosen such that Eq. (|35|) and Eq. ((26|) give estimation differing less than 5%. The canonical 
thermodynamical potential crff'J^j is then defined by Eq. (f22|) . with the identification Zp,^^ = Z^'^^^. 

Ensembles are equivalent if this function coincides with the entropy density (Tp^^^{p) defined by the Legendre 
transform of the grandcanonical partition sum Eq. p^ . More precisely, ensemble equivalence demands that the 
thermodynamic quantities as calculated from the canonical partition sum, 

1 acrf " 



dp 



-can 



dp 



(36) 



coincide with the corresponding grandcanonical quantities. The canonical thermodynamics, in the same thermody- 
namic conditions as in Fig. [TJ is displayed in Fig. [5] We can see that the canonical calculation allows to interpolate 
between the dense and diluted branches observed in the grandcanonical ensemble as expected. However the inter- 
polation is not linear, meaning that the chemical potential continuously varies as a function of the density. The 
discontinuity in the entropy slope at high density leads to a jump in the intensive observables close to the saturation 
density, in complete disagreement with the grandcanonical solution Even more interesting, the entropy presents 
a convex intruder, the behavior of the equations of state is not monotonous and a clear back-bending is observed, 
qualitatively similar to the phenomenon observed in first order phase transitions in finite systems 

The canonical partition sum Eq. ([50)) is only defined for finite values of the total number of particles A. The doubt 
therefore arises that this non-trivial behavior and the similarity with the finite systems thermodynamics might be 
due to a non achievement of the thermodynamic limit. Fig. [3] shows that this is not the case, and the canonical 
calculation is convergent. In this figure the chemical potential (a), pressure (b) and cluster size distribution (c) are 
represented for different values of the total number of particles A at a given density inside the in-equivalence region. 

We can see that indeed a very large size has to be used before this limit is achieved. As a rule of thumb, the 
total system size has to typically be approximately ten times bigger than the most probable cluster size in order 
to have convergent results, meaning that the equivalent of a Wigner-Seitz cell contains in average ss 10 dominant 
clusters. This can be understood considering that at finite temperature the distribution is very large, and the often 
used single-nucleus approximation [l^ is not realistic. 

As it can be seen in Fig. [3J the grandcanonical equilibrium prediction for this thermodynamic condition would 
correspond to a macroscopic liquid fraction in equilibrium with essentially free particles (dashed line). An explicit 
computation of the coexistence region in the canonical ensemble shows that this is not the case, as the cluster sizes 
do not scale with the total system size A. As one may notice, this liquid fraction is replaced by a finite, though large, 
nucleus with a characteristic radius of the order of only 5 — 10 fm. 

This finding is in agreement with all the theoretical microscopic modelizations which are naturally elaborated inside 
a single Wigner-Seitz cell within a fixed number of particles 3^- 1^ ■ AH these canonical studies agree in predicting 



that matter is clusterized for all subsaturation densities and the cluster size and composition evolve continuously 
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FIG. 2: Constrained entropy (a), pressure (b),(c) and chemical potential (b)(d) evaluated in the canonical ensemble in the 
same thermodynamic conditions as in Fig. 1. The pressure and chemical potential discontinuity is indicated by a dotted line. 



with the density, which is incompatible with a inodelization based on phase coexistence I40| where only the relative 
proportion of the two phases varies through the phase transition. 

As it has been argued in Refs. [l^, the quenching of the phase transition is due to the high electron inconipressibility. 
Because of the charge neutrality constraint over macroscopic distances, an intermediate density solution given by a 
linear combination of a high density homogeneous region and a low density clusterized region would imply an infinite 
repulsive interaction energy due to the electron density discontinuity at the (macroscopic) interface [16[ . 

This argument requires electrons to be completely incompressible. Since the electron incompressibility, while high, is 
not infinite, it could have been argued that a slightly modified coexistence region should be observed, where the liquid 
fraction would be constituted by large but still mesoscopic clusters, such that the interface energy would not diverge. 
The comparison between canonical and grandcanonical shown in Fig. [TJ [21 [21 demonstrates that this is not true: the 
presence of microscopic, instead than macroscopic, fluctuations, qualitatively change both the thermodynamics and 
the composition of matter. 
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FIG. 3: Convergence study of the canonical ensemble. The baryonic chemical potential (a), baryonic pressure (b) and cluster 
multiplicity per unit volume as a function of the cluster size are represented at a representative density inside the ensem- 
ble in-equivalence region, by varying the total system size. Thick line in panel (c): cluster multiplicity distribution in the 
grandcanonical ensemble. 
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FIG. 4: (a) Cluster distributions as a function of the density in the ensemble in-equivalence region; (b) mass fraction of 
unclusterized matter as a function of the total baryonic density, in the same thermodynamic conditions as in (a). The dashed 
line gives the saturation density of nuclear matter. 



III. PHENOMENOLOGICAL CONSEQUENCES OF ENSEMBLE IN-EQUIVALENCE 

It is important to remark that in the density region where the grancanonical ensemble is defined, the predictions of 
the ensembles coincide. The in-equivalence is observed in the intermediate density region, where the grandcanonical 
first order phase transition is not observed in the canonical ensemble. Since in this density domain the grandcanonical 
ensemble is not defined, there is no ambiguity on which ensemble should be chosen, meaning that the phenomenology 
of star matter has to be described with canonical thermodynamics. 

A closer look at the cluster distribution in the in-cquivalence region can be obtained from the left part of Fig. 
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in which displays the cluster distributions as a function of the density in the ensemble in-equivalence region. We 
can see that this distribution varies continuously, as it is expected physically, with very heavy nuclei and a very 
wide size distribution at the highest densities corresponding to the inner crust. These very massive nuclei disappear 
at a density close to saturation density, which defines the density corresponding to the crust-core transition in the 
model. This sudden process is at the origin of the non differentiable point observed in the entropy in Fig. [21 The 
resulting discontinuity in chemical potential and pressure is therefore a physical effect in the framework of this model. 



A word of caution is however necessary. It is well known from star matter literature j41| that close to saturation 
density deformed extended nuclei are energetically favored, the so called "pasta" phases. We expect that adding a 
deformation degree of freedom to the cluster energy functional would smooth this discontinuity. 

The right part of Fig. |4] represents the mass fraction of unclusterized matter (free nucleons and homogeneous 
dense matter) as a function of the density. Homogeneous matter dominates at the very low densities which physically 
correspond to the neutron star atmosphere, while clusters become increasingly dominant at higher density, until they 
melt into the homogeneous liquid core. Even in the absence of a first order phase transition, a very sharp behavior is 
obtained defining a relatively precise value for the crust-core transition density . 

Besides the relevance of the issue of ensemble in-equivalence from the statistical physics viewpoint, it is interesting 
to remark that the use of grandcanonical thermodynamics can lead to important qualitative and quantitative dis- 
crepancies in the computation of different physical quantities of interest for the astrophysical applications. This is 
demonstrated in Fig. [SJ which shows the cluster distribution for a chosen thermodynamic condition ( temperature 
T = 1.6 MeV, baryonic density p = 3.3 • 10^^ gcm~'^, proton fraction Yp = 0.41) which is typical for the dynamics 



of supernova matter after the bounce and before the propagation of the shock wave [Uj, [12] ■ We can see that the 



dominant cluster size is around A = 60, which is a particularly imp ortant size in the process of electron capture which 



determines the composition of the resulting neutron star [421 - 14 7] . It is clear that it is very important to correctly 
compute the abundances of such nuclei. 

Conversely in a grandcanonical formulation, as the widely used nuclear statistical equilibrium (NSE) [isl - Eil l , these 
partitions are simply not accessible as they fall in the phase transition region. An approach consisting in taking 
the mctastable grandcanonical prediction and considering only the nuclei of size such that the chosen total density 
is obtained, as in Ref. is shown by the dashed line in Fig. [5] It is clear that such approach completely 

misses the correct cluster distribution. Alternatively, hybrid canonical-grandcanonical formulations are routinely 
used in the astrophysical community [H, S, . Such approaches do not share the drawback of grandcanonical NSE 
calculations, but they always introduce artificial Maxwell constructions to fill the high density part of the equation of 
state. Moreover, they never address the fluctuations in the cluster composition, clusterized matter being modelized 
by a single representative nucleus. It is clear from Fig. [S| that this approximation is highly questionable at finite 
temperature, where distributions arc wide, the largest cluster docs not coincide with the average one nor with the 
most probable. The model we have presented overcomes all these problems. It is clear that many improvements 
are still necessary in this model before it can be considered as a reliable quantitative prediction for astrophysical 
simulations. Both improvements on the cluster energy functional and inclusion of deformation degrees of freedom 
are in progress to this aim. However we believe that the main results presented in this paper, namely the absence of 
phase transition due to Coulomb frustration and the dominance of microscopic clusters with a large and continuous 
distribution of size extending over most of the subsaturation region, are general results which will not change with a 
more sophisticated model. 



IV. CONCLUSIONS 



To conclude, in this paper we have shown that dense stellar matter as it can be found in core-collapse supernova 
and in the crust of neutron stars is a macroscopic physical example of ensemble in-cquivalcncc. A first-order phase 
transition is observed in the grandcanonical ensemble, but when the region corresponding to the discontinuity is 
explored explicitly constraining the density in the canonical ensemble, the macroscopic dishomogeneities associated 
to phase coexistence are seen to be replaced by microscopic dishomogeneities leading to cluster formation. As a 
consequence, the transition observed in the physical system is continuous. This phenomenon is due to the long 
range Coulombic interactions which quench the phase transition, thus giving rise to a thermodynamics qualitatively 
similar to the one of finite systems including thermodynamic anomalies. Specifically, the relation between density and 
chemical potential is non monotonous implying negative susceptibility. 

Accounting for this specificity can have sizeable effects in the computation of different quantities of interest for 
supernova dynamics. 

This paper has been partly supported by ANR under the project NEXEN and by IFIN-IN2P3 agreement nr. 07-44. 
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FIG. 5: Comparison between canonical (full line) and grandcanonical (dashed line) predictions for the cluster distribution in a 
specific thermodynamic condition relevant for supernova dynamics. 
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